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Abstract 

We continue the analysis of our previous articles which were devoted to type-I 
parametric down conversion, the extension to type-II being straightforward. 
We show that entanglement, in the Wigner representation, is just a corre- 
lation that involves both signals and vacuum fluctuations. An analysis of 
the detection process opens the way to a complete description of parametric 
down conversion in terms of pure Maxwell electromagnetic waves. 
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1 Introduction 



The theory of parametric down conversion (PDC) was treated, in the Wigner 
formahsm, in an earher series of articles|jl], 0, |^. There we showed that, 
provided one considers the zeropoint fluctuations of the vacuum to be real, 
the description of radiation is fully within Maxwell electromagnetic theory. 
Effectively, because the Wigner function maintains its positivity, we can say 
that quantization is just the addition of a zeropoint radiation, and there is 
no need for any further quantization of the light field. In the present article 
we show that the same result extends, without any difficulty, from the type-I 
PDC case to the type-II situation. 

There seems to be a widespread reluctance to accept the reality of the 
vacuum fluctuations, in spite of the fact that they appear, quite naturally, in 
the Wigner function of the vacuum state. We remark that such fluctuations 
have been taken seriously, within a certain school of thought, throughout 
the entire history of the quantum theory, following the formulation of Max 
Planck, originating in 1911[^. Of course, it is true that, integrated over 
all frequencies, they give us a vacuum with infinite energy density; why 
then are all photographic plates not blackened instantaneously? But all 
photodetectors, including even our own eyes, are very selective, not only as 
regards the frequency, but even also the wave vectors, of the light components 
they analyze. This is especially the case with the detectors commonly used 
in PDC experiments. So, there is a noise to subtract, but it is not infinite! 

In our previous articles we indicated how the noise subtraction is made, 
according to the Wigner formalism, and showed how this subtraction is re- 
lated to the standard calculating device, of normal ordering, used in the 
Hilbert-space formalism. Here we extend this analysis, in an informal man- 
ner, showing that, if we take into account the fact that all detectors integrate 
the light intensity over a large time window, the process of light detection, 
like that of light propagation, may also be described entirely in terms of 
real waves and positive probabilities. We are then able to see that, in terms 
of a purely wave description, the highly problematic concept of "entangled- 
photon" states of the field loses all its mystery. Entangled photons are just 
correlated waves! The only reason this description has taken so long to ma- 
ture is that the word "classical", in reference to the light field, is restricted 
in its application to Glauber-classical states 0. A discussion of the differ- 
ence between classical and nonclassical effects has been given in Ref. . The 
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states which are produced when a nonhnear crystal interacts with a coher- 
ent incoming beam, and, of course, simultaneously with the vacuum, may 
be described using classical Maxwell theory, but there is correlation of the 
outgoing light beams both above and below the zeropoint level. 

2 General description of parametric down con- 
version in the Wigner representation 

Type-I parametric down-conversion, in which the correlated signal and idler 
beams have the same polarization, has been recently studied within the 
framework of the Wigner function |^, ^ ^. The formalism is almost identical 
in the case of type-II PDC, in which the correlated beams leaving the non- 
linear crystal are orthogonally polarized. The process of type-II PDC can be 
formalized in analogy with the classical Hamiltonian of [jl| 

j=o,e k 

+ {iWV /(k, ^')exp{-iujpt) alkali,, + c.c), (1) 

k,k' 

o(e) refers to the ordinary (extraordinary) rays. We have taken the origin 
of the coordinate system at the center of the crystal, and treated the pump 
beam as an intense monochromatic plane wave represented by 

V(r, t) = {V{t)exp[i{kp ■ r - Ljpt)] + c.c.) u, (2) 

u being a unit vector perpendicular to kp. As the coherence time of the laser 
is large in comparison with most of the times involved in the process, we shall 
consider V{t) as a constant, g' is a real coupling constant defined so that 
the product g'V has dimensions of frequency, and /(k, k') is a dimensionless 
symmetrical function of the wave vectors inside the crystal. This function, 
which is related to the function h{]<., k') introduced in equation (8) of refer- 
ence 1^ , is different from zero only when the following matching condition is 
fulfilled 

kp ^ k + k'. (3) 

As is well known 0, there is in addition a matching condition for fre- 
quency that is fulfilled much more rigorously, namely. 
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ujp — uj\^ + a;k'- 



(4) 



On the other hand, ctok (ctek') is the field amphtude for the mode with 
wave number k (k') corresponding to the ordinary (extraordinary) field, 
which is represented as a sum of two mutually conjugate complex c-numbers 

E,(r,t) = E(+)(r,t)+E(-)(r,t), (5) 

Ej^V.^) =^Y^ (^^) eikaik(t)exp(ik ■ r), j = o, e, (6) 

where is the normalization volume and e^k is a polarization vector. Eqs. 
(ID and (H) correspond to the Heisenberg picture, where all time dependence 
goes in the field amplitudes Oi*]^{t) and aj\^{t). For a free field this dependence 
has the form 

ajk(^) = ajk(0)exp(-iu;jk^), (7) 

but for interacting fields it is complicated and contains all the dynamics of 
the process. 

The evolution of the Wigner field amplitudes aj\^{t) is directly given by the 
Hamilton (canonical) equations of motion taking \/haj\^{t) as coordinates and 
V^oc*j\Sf) as canonical momenta. For instance, we get for the extraordinary 
field amplitude Oek: 

dek = -^t^ekttek + q'V ^ /(k, k')exp(-zu;pt)a*k', (8) 

k' 

and a similar expression holds for the ordinary field amplitude by exchanging 
index e with o. 

In order to calculate aek('^) for all t we shall take into account that the 
amplitude aekit) evolves as a free-field mode before entering the crystal and 
after coming out. We shall integrate (§) from t = —At to t = 0, where At 
is the time taken for the radiation to cross the crystal. The initial condition 
is aek(— At) = a^^^\—At), where a^^'^\—At) is the field amplitude of the 
mode k in the incoming vacuum field. 

To second order in the coupling constant g', that is taking the second 
term of the right side of (|^) as a perturbation and retaining terms up to 
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order g'"^, we get (putting g' /S.t = g) 

k' ^ 
VI^P E E /(k, k')r (k', k")«[^(u;ok' + o^ek" - 0.,)] 

k' k" ^ 

rAt, (vac). 



xn[— (^ek"-^ek)]«^;^:^^(0) ; (9) 



where 



u{x) = exp{ix), (10) 

oc 

Equation (|^) implies k" k in the second order contribution to (|^). 
In the derivation of {^) we have taken into account that 

a^J^''\o) = a^J^''\—At)exp{—iuJe]iAt), 

«1™'^(0) = a:t''\~At)expituj,^At). (11) 
After t = 0, aek{t) evolves free-field mode 

Oiekit) = aek(0)exp(-it^ekt). (12) 

Now, let us consider two narrow correlated beams called "ordinary" and 
"extraordinary", with average frequencies Uq, uje, and wave vectors ko, kg 
respectively, fulfilling the matching conditions 

uJo + uJe = ujp ; ko + kg = kp. (13) 

Both light beams contain frequencies within a range between ujjram and cu^max 
(j = o, e), wave vectors whose transverse components are limited by a small 
upper value, and orthogonal polarization vectors which are practically inde- 
pendent of the wave vectors, that is 

^ ^ I XA'^ I ^ ^imin 

eek = Ce , eok = Co ; ee ■ eo = 0. (14) 
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We also substitute slowly varying fields F(+)(r,t) (F(-)(r,t)) for the ampli- 
tudes E(+)(r,t) (E'^~)(r, t)), the relation between them being 



f1+)( 



E 



ttk(0)exp(zk ■ r)exp[i{Ljj — ujjkjt] 



ej (j = o,e), (15) 



where Uj is some appropriately chosen average frequency midway between 
tUmin and tUinax (see(p!4D). The square brackets in the summation symbol 
indicates that the sum is restricted to the set of k pertaining to the j^beam. 
It is easy to obtain the amplitude ¥j^\r B,t) in terms of the amplitude 

F^j^\rA,t) at another point of the light beam |l|]. We find 



rB,t) = F'j^\rA,t- '-^)exp{iujj'-^) {j = e,o) 



Tab 



rABs 



(16) 



where tab = rs - and tab = \t^ab\- 

From expressions (0), and (|T5|) we obtain 



([1 + ^2|\/|2j]FW(™^)(r,t) + ^yGFi-)(™^)(r,t)) 



(17) 



F(+)(r,t) = FW(r,t)e. 
'[1 + /|r|V]FW(™^)(r,t) + ^rGFj-)(™^)(r,t)) e^. 

Here F^™'^) and F^™'^) are the incoming vacuum fields and Fg (Fq) the out 
going extraordinary (ordinary) fields — see Fig. 1. We have 



(18) 



Fi+)(-^)(r,t)=^5: 



ek 



2L3 



exp{ik ■ r)exp[{uje - uJe]i)t]a''J^''\o) , (19) 



and similarly for _Fj+)(™'=). G and J are linear operators which are defined 
as: 



GFi-)(-^)(r,t)=.E(|^^ 



exp{ik ■ r)exp[i{uje - u;ek)^]/5k, 



(20) 
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with 

/3k = E /(k, k')«[^K - ^ek - uj,^>)]atr\0), (21) 



and 



j^(+){vac) = , ^ " expizk ■ r)exp[z{u, - c.,k)t]7k, (22) 



with 



7k = E E /(k,k')r(k',k>[^(a;,k' + ^ek"-o;p)]x 

[k']c [k"]e 

X.[^(.ek"--ek)]«ik"^(0). (23) 

From (^) we see that the outgoing extraordinary beam, to order g'^, con- 
sists of three parts: i) A zeropoint radiation with amphtude which 
passes through the crystal without any change; ii) a radiation produced by 
the nonhnear interaction (mediated by the crystal) between the laser beam, 
with amplitude V, and the zeropoint radiation, with amplitude Fj;~^^'^^^\ en- 
tering the crystal in the direction of the ordinary beam; iii) one part which 
just modifies a little (to order g"^) the amplitude F^~^^'^'^^'^\ The ordinary beam 
is constituted in a similar manner. 

Now, let us consider the correlation properties of the fields, which are 
identical to the ones calculated in a previous work 0: 

a) Autocorrelations. 

Taking the extraordinary field F^+)(r,t) = Fj+^(r,t)ee at a point r and 
times t and t', we have: 

{F^+\r,t)Ft\r,t')) - (Fi+)(™^)(r, t)Fi-)(-^)(r, t')) = 

2/|y|2(GFi-)(™^)(r,t)G*FW(™=)(r,t')) = g^\V\^i^e{t' - t), 

{Ft\r,t)F^^\v,t'))=0. (24) 

Here "()" means an average using the Wigner function in the vacuum state 
as probability density. /ie(t — t') is a correlation function which goes to zero 
when |t' — t| is greater than the correlation time of the extraordinary beam, 
Tg. Similar expressions hold for the ordinary field by exchanging the indices 
"e" and "o". 
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b) Crosscorrelations. 

Taking the extraordinary (F^^^ (r, t) = Fj;~^^ (r, t)ee) and ordinary (F^+) (r, t) = 
Fj+^(r, t)eo) fields at the center of the crystal r = r' = and times t and t', 
we have: 

{F^^\0,t)F^+\0,t')) = gVuit'-t). 

{F^^\0,t)F^-\0,t')) = (Fi-)(0,t)FW(0,t')) = 0. (25) 

Here z/(t' — t) is a function which vanishes when \t' — t\ is greater than the 
coherence time between the extraordinary and ordinary beams. From (|25|) it 
is possible to derive all crosscorrelations at different points r 7^ r' by using 

Finally, the quantum theory of detection in the Wigner representation 
gives us the following results for single and joint detection probabilities: 

a) Single probability. 

The following result is a general expression for calculating single proba- 
bilities per unit time in the Wigner representation: 

Pi(ri,t)(x(/(ri,t)-/o(ri)), (26) 

where /(ri,t) = |E(+)(ri,t)p = |F(+)(ri, t)p, and /o(ri) is the intensity of 
the vacuum field at the position of the detector. 

b) Joint probability. 
It can be proved that in PDC experiments 

Pi2(ri, t; V2, t + r) oc ({/(ri, t) - /o(ri)}{/(r2, t + r) - /o(r2)}). (27) 

By taking into account that the Wigner fields amplitudes are Gaussian pro- 
cesses, and neglecting fourth order terms in g, we have |Q: 

P,,(ri, t; r2, t + r) cx 5: 5: | (F^ (n, t)^^ (r^, t + r))\\ (28) 

A A' 

where A and A' are polarization indices. 
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3 Tests of Bell's inequalities using polariza- 
tion correlation 



Most experiments to test Bell's inequalities with nonlinear crystals performed 
hitherto have used type-I parametric down-conversion in which the two cor- 
related beams have the same polarization. In |1|, 0, experiments of this 
kind were analyzed in the Wigner function formalism. However, more re- 
cent experiments, using type-II phase matching, provide a more direct way 
to generate "entangled-photon" states. Type-II experiments are themselves 
of two types. In the first, that is coUinear type-II PDC, the crystal is ori- 
ented so that the ordinary and extraordinary radiation cones are mutually 
tangent in the direction of the pumping beam. To date, nearly all type-II 
experiments have used coUinear phase matching [|]. On the other hand [^j, 
in noncoUinear type-II phase matching, the two cones intersect along two 
directions, and this gives rise to an entangled state in the polarization (see 
Fig. 2). It has been claimed that such a source produces true entangled states, 
capable of violating Bell's inequalities. 

The experimental outline is shown in Fig. 3. The two beams "1" and "2", 
in which the ordinary and extraordinary cones intersect, are selected and sent 
to two polarizers Pi and P2 oriented at angles 0i and 02 with respect to the 
polarization of the extraordinary ray. Coincidence rates were measured as 
functions of angles 0i and 02. In additional optical devices, that is half- 
and quarter-wave plates, were used in order to produce four different Bell 
states, but we shall confine our analysis to just one of these states, namely 
the one which uses no additional devices. 

Let us see how the entangled state is represented in the Wigner formalism. 
The two beams, coming out of the crystal along the directions where the 
ordinary and extraordinary cones intersect, are given by 

FW(0,t) = FW(0,t)i + Fi+)(0,t)j, 

F^^\0,t) = Fi+)(0,t)i' + FW(0,t)j', (29) 

where i, i' represent the polarizations of the extraordinary beams and j, j' 
the polarizations of the ordinary beams. The essential point is that the 
extraordinary component, F^^\ of the first ray and the ordinary component, 
of the second ray are conjugated, and therefore correlated. Similarly, 
and are correlated, but F^+'> (Fj+)) is uncorrelated to (F^^^). 
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When a polarizer oriented at angle 0i to the horizontal is placed in front 
of the detector Di, the field at Di (placed at ri) at time t is 

FW(ri,t) = [Fi+)(ri,t) ■ (cos0ii + sin0ij)](cos(/)ii + sin^J) 

= exp[iuj{d/c)][F^^\0,t- ^) cos(/)i + Fi+^(0,t - ^) sin^i] 

(cos 0ii + sin 0ij) (30) 

Here we recall that the action of a polarizer is to project the electric field 
vector on the polarization direction 

In the same way, we write for the field at the detector D2 (placed at r2) 
at time t + t 

F(+) (ra, t + T) = [F^+^ (ra, t + r) ■ (cos ^ai' + sin 02j')] (cos ^21 + sin 02j') 

d d 
= exp[iu{d/c)][F^;^\0,t + r - -) cos 02 + F^+\0,t + r - -) sin^a] 

(cos 02 i' + sin 02j') (31) 



In order to calculate the joint probability we combine Eqs. (|28|) , (|30|) , 
and (0), and take into account the correlation properties of the fields given 



by (^^ and (|^). After some easy algebra and an integration of Pwi^) over 
the detection window, we obtain the coincidence probability 

Pi2 = Ksin2(02 + 0i), (32) 

K being a constant. This expression is similar to the one obtained in 0, and 
corresponds to 100% contrast. This type of correlation is usually claimed to 
violate a Bell inequality (but see the discussion section). In the actual exper- 
iment a violation of the inequality by 100 standard deviations is reported. 



4 The detection problem in the Wigner rep- 
resentation 

The detection probability in quantum optics is usually written in terms of 
the normally ordered expression 
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/•At „ 

Poc / {E^-\t)E^+\t))dt, (33) 
Jo 

where E^'^^ (t) is the Heisenberg operator of the electric field at the detector, 
E^~\t) its hermitian conjugate and At is the detection time window. For 
simplicity we consider a point detector and, therefore, ignore the standard 
volume integral in (p^). 

When we pass to the Wigner representation, the normally ordered expres- 
sion should be written in terms of a symmetrically ordered expression minus 
a commutator. Then we may replace (see Eg. (p6|) ) the Heisenberg operators 
by random wave amplitudes and, after some rather trivial algebra, we get 

i-At 

P (X / {I{t) - Io)dt, (34) 





where I{t) = \E{t)\'^ is the intensity of the field arriving at the detector and 
/o is a constant corresponding to the average intensity of the zeropoint, i.e. 
the intensity that would arrive at the detector if all light sources, such as 
lasers, were switched off. 

If the Wigner function is positive definite, I{t) may be interpreted as 
a stochastic process, which makes possible a wavelike interpretation of the 
propagation of light. This is the case in all experiments involving parametric 
down conversion. However, there remains a problem for a wavelike inter- 
pretation of the detection process, because I(t) — Iq is not positive-definite. 
This means we cannot assume that I{t) — Jq is proportional to a detection 
probability. Nevertheless, it is easy to show that the average {I(t) — Iq) is 
nonnegative-definite. Now we consider the usual case where I{t) is a sta- 



tionary stochastic process. It then follows (see Ref.|]TO[, page 49) that it is 
ergodic. We could consider substituting time averages for ensemble averages, 
thereby obtaining 

P oc {lit) - /o) = hm i /^[/(t) - Io\dt. (35) 

T->oo 1 Jq 

This shows that, for every member of the ensemble of stochastic wave am- 
phtudes (except for a subensemble of zero measure), the time-averaged pho- 
todetection probability precisely equals the ensemble averaged one (/(t) — Jq). 

In practice we do not have a time-average but something which is almost 
equivalent, namely an integration over the detection window At. A typical 
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detection window lasts more than one nanosecond, so that the dimensionless 
quantity tuAt, where uj is the frequency of visible light, is of the order of 10'^. 
It is true that taking the limit T — oo is not the same as taking T 10'' (in 
dimensionless units), but the difference must be small. Since the right hand 
side of ( p5|) is nonnegative-definite, this means that the finite-time average 

—j^ m-i,]dt (36) 

takes negative values only with a very small probability. Now let us modify 



the standard detection probability of Eq.(34) so that 




P cx / [/(t) - h]dt , (37) 



where the notation {}_|_ indicates that we replace the contents of the brack- 
ets by zero if their value is negative. Then ( PTf ) will predict rather more 
photocounts than the standard quantum detection theory, but it seems not 
unreasonable to assume that these additional counts correspond to a part of 
the dark rate at the detector. As a matter of fact, the quantum theory of 
detection, leading to (^3|), involves first order perturbation theory. Therefore 
we may assume that quantum theory also predicts some dark background in 
photocounters when higher-order processes are taken into account, because 
the detector may be activated by vacuum fluctuations. A detailed study of 
such a background would have to enter fully into the electronic band struc- 
ture of the photodetector material |11[]; a treatment based on the photoelectric 



effect for single atoms is clearly inadequate. 



5 Discussion 

In the preceding sections we have outlined a theory of both the propagation 
and detection of light which is consistently local realist in the sense defined, 
for example, by Clauser and Shimony|jl2|. In fact, eq . (^Tf ) has the standard 
form introduced by Bell in his definition of local hidden variables models, in 
particular 

Pi2 = (Pi(A,0i)P2(A,02)) = / p(A)Pi(A,0i)P2(A,02)c^A, (38) 
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where the amphtudes a^^k involved in /(r,t), via E{r,t) (see eq.(||) and the 
hne after eg . (|2^) ) play the role of the hidden variables A. On the other 
hand, our theory is also in almost perfect one-to-one correspondence with 
the standard Hilbert-space theory, the only difference being the modification 
in the detection probability which we have proposed in Eg . (|37|) . Apart from 
a small dark rate, which is, in any case, a feature of all real experimental 
situations, this theory gives singles and coincidence rates agreeing with the 
standard theory . Therefore we have here an apparent contradiction. On 
the one hand, a Bell inequality is violated according to Ref.0, which implies 
that no local hidden variables model exists for the experiment. On the other 
hand we have an explicit local hidden variables model, namely the guantum 
model in the Wigner representation with the modification of (|T7D. What is 
the resolution of this apparent contradiction? The reason is that the Bell 
ineguality which is actually violated in the experiment is not a genuine Bell 
ineguality derived from the assumptions of realism and locality alone (like 
ineguality (4) of Clauser and Horne flSl), but a homogeneous ineguality 
involving additional assumptions (like ineguality (11) of Clauser and Horne). 
One of these assumptions, rather than local realism, is what is violated in 
our LHV model. 



The present authors have been insisting for many years now ||Tj] that 
these auxiliary assumptions are not only unreasonably restrictive but also 
incorrect. We believe that the results reported in the present article vindicate 
our point of view. 

As is well known a genuine Bell ineguality cannot be tested in experiments 
with visible light, due to the low efficiency of the detectors presently available. 
The conventional wisdom is that this is just a technical problem that will 
be solved in the near future. However, the existence of an LHV model for 
the guoted experiment which does not rest upon the low efficiency of the 
detectors, but on the existence of some unavoidable amount of dark rate (see 
Eg.(P7|)) shows that any future rehab le test of LHV theories should involve 
detectors having both high efficiency and low dark rate. 

Now, we turn to the interpretation of entanglement in the Wigner repre- 
sentation. We saw (see Egs.(p5D and (^)) that "photon entanglement" is just 
correlation between two light beams. Then, what is the difference between 
"classical" correlation and entanglement? In guantum optics what is usually 
called "classical" is light having a (Glauber) P-representation which is posi- 
tive definite, and "classical correlation" usually means a correlation between 
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the P-distribution functions of the two hght beams. Here we see that "local 
realist theories" is a class much bigger than standard "classical theories" . In 
particular, the quantum model following from the Wigner representation in- 
terpreted as a (positive) probability distribution allows for correlations much 
stronger than "classical correlations" . As in our functions F^~^^ and Fj;^^ of 
Eq. (p5|) they involve correlation of the zeropoint part of the electromagnetic 
field, and this is "entanglement". 

Of course, the obvious objection to these arguments is that they cannot 
be extended to cases where the Wigner function is negative. Elsewhere we 
have conjectured that this never happens in actual experiments. We refer to 
our publication p for details. 
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Figure captions 

Fig. 1. The process of parametric down conversion. 

Fig. 2. Polarization entanglement in noncoUinear type-II down conver- 
sion. 

Fig. 3. Tests of Bell's inequalities using noncoUinear type-II down 
conversion. 

Figures have not yet been included, as the authors have to master the 
appropriate technology first. We hope it will not take very long! 
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